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Abstract 

A superconductor with 4-fermion attraction, considered by Mackowiak and Tarasewicz is modi- 
fied by adding to the Hamiltonian a long-range magnetic interaction V between conduction fermions 
and localized distinguishable spin 1/2 magnetic impurities. V has the form of a reduced s-d in- 
teraction. An upper and lower bound to the system's free energy density f(H,f3) is derived and 
the two bounds are shown to coalesce in the thermodynamic limit. The resulting mean-field equa- 
tions for the gap A and a parameter y, characterizing the impurity subsystem are solved and the 
solution minimizing / is found for various values of magnetic coupling constant g and impurity 
concentration. The phase diagrams of the system are depicted with five distinct phases: the nor- 
mal phase, unperturbed superconducting phase, perturbed superconducting phase with nonzero 
gap in the excitation spectrum, perturbed gapless superconducting phase and impurity phase with 
completely suppressed superconductivity. 

PACS numbers: 74.20.Fg, 74.25.Bt, 74.25.Dw 
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I. INTRODUCTION 



The properties of superconductors can be drastically changed by adding new elements, 
e.g. the superconducting transition temperature T c decreases in most superconductors, when 
magnetic impurities are added 1 . On the contrary, Bednorz and Miiller discovered that copper 
oxide doped with lanthanum and barium exhibits high T c pa 35 K-. 

The effect of magnetic impurities on superconductors was studied by Abrikosov and 
Gor'kov^. Using Green's functions they not only explained the strong decrease in T c in the 
presence of magnetic impurities, but also predicted "gapless behaviour" of superconductors 
and complete destruction of superconductivity above critical concentration of impurities. 

The problem of lowered T c in the presence of magnetic impurities was also discussed 
by other theorists, e.g. Nakamura^ and Suhl et al.- explained this effect by treating the 
s-d interaction V s _^r as an additive term in the total Hamiltonian, which perturbes a BCS 
superconductor-. Balseiro et al.- studied a BCS superconductor perburbed by magnetic 
impurities interacting via a nearest neighbour Heisenberg potential. The resulting phase 
diagrams comply qualitatively with experiment. Simon and others- analyzed the disappear- 
ance of the Kondo effect in d-wave superconductors. Openov^ extensively discussed the 
question of reduced specific heat jump at T c in superconductors containing magnetic and 
nonmagnetic impurities. 

Our objective here is to study, the effect of magnetic impurities on a superconductor in 
a volume |A|, with 4-fermion attraction of the form 



and a k+J a k _ are fermion annihilation operators, whereas G kk , is real, symmetric, invariant 
under the time reversal k — > —k or kS — > — and nonvanishing only in a thin band close to 
the Fermi surface, viz., 




(1) 



where 



frfc _ a k+ a k-i 



G 



kk' ~ 



G x (k)x(k'), 



G > 0, 



where x(k) denotes the characteristic function of the set 



S = {k: // - 5 < e k < /i + 5} 



h 2 k 2 



2m 
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A superconductor with an interaction of this type was studied by Mackowiak et 
al.— li^iMiM. Tarasewicz et al.— demonstrated that 4-fermion attraction appears in sixth or- 
der expansion terms of Frohlich's transformation^. So far, there is not much experimental 
evidence, which could confirm existence of fermion quartets generated by V^f in supercon- 
ductors. An exception is the discovery of half-/i/2e magnetic flux quanta in SQUIDs 17 . Flux 
quanta with values close to h ~ 6/e and h ~ 8/e have also been observed. Accordingly, we 
can expect electrons to bind into quartets, sextet's, octets, and so on. 

In this preliminary step of the investigation of a superconductor with 4-fermion attraction 
perturbed by magnetic impurities, we assume the perturbation implemented by the localized 
distinguishable spin 1/2 magnetic impurities to be a reduced long-range s-d interaction, 
which involves only the z-components of the impurity and fermion spin operators. The 
reason for this simplification is that the thermodynamics of the resulting Hamiltonian H = 
Hq + V±f + V admits a mean- field solution which improves with decreasing impurity density. 
Furthermore, this solution is thermodynamically equivalent to the one obtained for H with 
a Heisenberg type reduced s-d interaction 

f 
N 



ka 

°a± = °~ax i l(T ayi 

replacing V (Ref.Q, Sec. 6.2.5). (Similarly, the thermodynamics of classical superconductor 
can be explained in terms of a reduced BCS interaction, whereas a gauge-invariant theory 
of the Meissner effect requires a more general pairing potential.) The reduced form of Vh, 
obtained by rejecting the sum V^y in the s-d interaction V^_d, is an approximation which 

k^k' 

resorts to the fact that spin-exchange processes, and not momentum exchange processes, are 
primarily responsible for the Kondo effect caused by K-<r^- 
The interaction V has the form 
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V = -Tr^K- -n k+ )<r a , v a =\ , (3) 

denoting the number of magnetic impurities. The interaction which induces the super- 
conducting transition was chosen to be V4/. The resulting theory is relatively simple and 
allows to construct the system's phase diagrams. To this end, the system's free energy 
density f(H,/3) is bounded from above and below by mean-field type bounds, which are 



shown to be almost equal if the impurity density d and width of conduction band S are 
sufficiently small. The proof exploits methods developed by Czerwonkc 20 and Tindemans et 
al.— The resulting mean-field equations, for the gap A and a parameter y characterizing 
the impurity subsystem, are solved and the solution, which minimizes / is found for various 
values of magnetic coupling constant g and density d. These equations have different forms 
in the limit as temperature tends to OK, which depend on g. For sufficiently small values 
of g, which preserve the weak perturbation by magnetic impurities, the limiting form of the 
gap equation, as T — > 0, is the same as in the BCS theory. In the present paper the weakness 
of reduced s-d coupling is assumed. 

Phase diagrams of the system are depicted with five distinct phases: the normal phase, 
unperturbed superconducting phase, perturbed superconducting phase with nonzero gap in 
the excitation spectrum, perturbed gapless superconducting phase and impurity phase with 
completely suppressed superconductivity. An application of the theory to a superconductor 
with the electron-electron attraction of the form Vbcs + V4/, is under construction. 

II. THE HAMILTON! AN AND LOWER BOUND TO THE FREE ENERGY 

The full Hamiltonian of the system is 



where H = J2k<j £k n k<n wrfc h £fc = e k~^-> n ka = a la a ka, * s the kinetic energy of free fermions. 
In order to evaluate the free energy of the system 



it is convienient to separate the fermion and impurity operators. This can be done by 
exploiting the identity 



H = H' S + V, 



H' s = H + V 4 {, 



(4) 



T = -p-HnTrexp[-pH) 




ka 



ka 




ka 



±N\S\I, 



where I stands for the operator 



h ® • • • ® In, 
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and \S\ denotes the number of points in S. Let us now define 

H S = H + V 4i + ^N' 1 ^^- ~ 2%_% + + n k+ ) 



ka 



K + \g 2 \S\I, 



where 



Ho = H + V if + \g 2 N~ l J2( n k- ~ 2n k _n k+ + n k+ ). 
Then H can be written as 



ka 



H = H S + V, 



where 



9 c 

V ' = -^J2( n k-- n k + + °aY 



ka 



To linearize the interaction (J7J) we use the identity 

oo 

exp(a 2 ) = / exp(— |x 2 + V2ax)dx. 

v27r J 

— oo 

Then, exploiting the commutation relations 

[H s , n k _ - n k+ ] = [H s , a a ] = 0, 

the statistical sum of the system 

Z = Trexp(-/3#) 



Tr exp 



ka 



can be expressed as 



Z = Tr J Yl (dx ka exp \-\N(3x\ a + (3gx ka (n k - n k+ + a a ) 2 ) exp [ 

ka 

t\to N \ s \ r 



Z ka eX P 



\fiNx\ a + G ka {x ka ) 



where H s = Y Jk H sk = N 1 Y,ka H ska and 



Gkai^ka) = hiTrexp(-/3iV 1 H ska + Pgx ka (n k __ - n k+ + a a 



Following Pearce et al.— , let us now add and subtract \N f3(,~ l x\ a1 with £ > 1, in the 
exponent of the integrand on the r.h.s. of Eq. (fill) , and bound the resulting expression from 
above: 

N\S\ 



Z <Y[ maxexp [-fA^C^L + G 



ka 

x / Qexp 

k'a' 



ka \ X ka 



2tt 



(12) 



iN{3xl w (1 - C 1 ) 



dx h , ,. 



Eq. ffl2l now yields a lower bound to the system's free energy \A\f(H,/3) = —(3 1 lnZ: 

lAI/^^minffC 1 ^ (13) 



The necessary condition for the minimum on the r.h.s. of Eq. (1131) is 

Tr K- - n k+ + a a ) exp [~f3N~ l H ska + Pgx ka (n k _ - n k+ + a a )] 



C Nx ka = g 



Trexpf-^iV- 1 ^ + Pgx ka (n k _ - n k+ 



+ *«)] 



(14) 



for k e S, a = 1, . . . , N. 



III. UPPER BOUND TO THE FREE ENERGY 



An upper bound to the free energy of the system can be expressed in terms of the 
Hamiltonian h: 



h = H s + lNj2xl a -gJ2 X ka ( n k~ - n k+ + °a): 



ka 



where x ka € M, is a solution of the equation 

Nx ka = g {n k _ - n k+ + a a ) h 

with 



(n k _ - n k+ + a a ) 



Tr (n k _ - n k+ + a a ) exp(-{3h) 



k+ ' " a/h Tr exp(-/3h) 

The free energy of the system described by the Hamiltonian (Tl5|) has the form 



|A|/(M) = liV^^L-rMnTrexp 



ka 



~(3H S + Pg^, X ka ( n k- ~ n k+ + °a) 



ka 



lA^xL-ZTMnflTrexp 



ka 



ka 



(15) 



(16) 



(17) 



-PN^H^ + fax^ 



riu — n 



k+ 



(18) 



ka 



ka 
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Bogolyubov's inequality 

F(H,f])<F(h,P) + (H-h) h , (19) 
yields the relevant upper bound. Exploiting Eq. (TTBT) one obtains 

(H-h) h = (V) h -i N j24«+ N Yl x l = ( v ') h + 1 2 nJ2 4*> 

ka ka ka 

where 

(V') h = -WN- 1 £ ((K- " n k + Y) h + 2 <K- " % + ) °a) h + K) h ) • (20) 

ka 

The inequality Tr(pA 2 ) > (Tr(pA)) 2 , valid for any bounded self-adjoint operator A and 
density matrix p, implies 

(V) h < -hg'N- 1 ( n k- - H + + = -InJ2 x I- 

ka ka 

The thermal average (H — h) h is therefore bounded from above by and from Eq. ffl9l we 

get 

F(H,/3)<F{h,/3). (21) 
IV. COMPARISON OF THE UPPER AND LOWER BOUND TO f(H,/3) 



The form of Eqs. (fl4l) . (fl6l) shows that their solutions Xka do not depend on k, a. 
Furthermore, for |<^ _1 — 1| <C 1 both equations have almost identical solutions. Henceforth, 
they will be denoted by x. It follows now from the inequalities ffl3|) . f|2T]) and Eq. ffl8l) that 
the equality 

f(H,/3) = wmf(h(x),P), (22) 

X 

where 

h(x) = H s + \N 2 \S\x 2 -gx^ ( n k- ~ n k+ + °a) , 



ka 



holds up to negligible terms if the density of magnetic impurities and the width S of the 
conduction band S are sufficiently small. On these grounds we shall assume that the ther- 
modynamics of the original system is equivalent, under these restrictions, to that of h(x) 
where x is the minimizing solution of Eq. ffTol) . 



Refs. 



The thermodynamics of a system with a Hamiltonian of the form h(x) was studied in 

inn • 



In the next section this method is exploited to evaluate limmin f(h(x), (3) 



as |A| — > oo. 
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V. THERMODYNAMIC EQUIVALENCE OF MEAN-FIELD THEORY FOR h(x) 



The thermodynamic perturbation method of Bogolyubov et alM^ , extended by Czer- 



wonko in Ref. 



20, shows that in genera 
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lim f(h(x),/3)= Urn mm f(h (x), (3), (23) 

A|^oo |A|— >oo {h } 

where 

W=E^) + ^) + C„ (24) 

fc>0 

and 

h k( x ) = (4 + \g 2 ) Y ("far + n -ko) - 2A k{Pk + P*k) - 9 2 { n k- n k+ + n -k- n -k+) + U k( X ) + C kJ 

(25) 

U k (x) = -gxN(n k _-n k+ + n_ k _-n_ k+ ), C k = 2A k T k , C = ±\S\ (g 2 + N 2 x 2 ), (26) 

A k = l A l _1 Y G kk! T k>i h A x ) = Y h <*( x )> h ^ x ) = ~9 x \S\cr a , (27) 

y a 

with £ k = e k - //, l3 k = b_ k b k . 

The structure of U k (x) and h a (x) shows that h k (x) favours opposite alignment of impurity 
spins and those of conduction fermions. The interaction V therefore acts as a quartet- 
breaker. 

The difference 

Ah = h{x) - h (x) = -\A\- 1 Y G kk >B k Bl„ 

kkf 

where B k = f3 k — T k ,r k G Mr, yields a negligible contribution to the limit on the r.h.s. of 
Eq. (I2"3"j) . provided the constants r k are adjusted to satisfy the condition 

. . Tr (B k exp[-f3h k {x)]) , 
(B k ) = — V ~ r -\ = ]J = 0. 28 



According to Eqs. ( l23l) . ( |24l) the system, can be equivalently described in terms of hJx) 
and to evaluate f(h Q (x),[3) it suffices to diagonalize each h k . This was done in Ref. [13| by 
the method of Czerwonkc 20 . We briefly recapitulate the results. The operator h k acts in the 
16-dimensional space M k spanned by the vectors 

Kn 2 n 3 n 4 ) = {a* k+ ) ni (4_)™ 2 (a^)" 3 (a V)" 4 |0) , 



where = 0, 1; i = 1, 2, 3, 4. Let us now define the spin operator 

2S k = E E ar w ( 29 ) 

(p=±l cr=±l 

and seniorities 

A fc,a = n k,a - U -k,a, ° = ±> ( 30 ) 

and A fc(T commute with /i fc : 

[^,2Sj = [h k ,A ka ]=0. (31) 

The commutation relations f[3"Tj) enable diagonalization in the invariant subspaces of 
with fixed eigenvalues of the operators 2S k , A k+ , A k _ and h k . The space M k splits into nine 
such invariant subspaces, viz., 

— four 1-dimensional subspaces, which are spanned by the following eigenvectors with 
the corresponding eigenvalues 2s, A CT , E k of the operators 25^, A k a and h k — C k : 

|1010) 2s = 2 A + = l A_ = -l E k = 2£ k + g 2 + 2gNx 

|0101) 2s = -2 A + = -1 A_ = l E k = 2£ k + g 2 -2gNx 

|1100) 2s = A + = 1 A_ = 1 E k = 2^ k 

|0011) 2s = A + = -1 A_ = -1 E k = 24 
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four 2-dimensional subspaces, where the eigenvectors and eigenvalues are: 

|1000) 2s = 1 A + = l A_ = E k = i k +\g 2 + gNx 



0001) 2s = -1 A + = -1 A_ = E k = £ k + \g 2 - gNx 



0010) 2s = 1 A + = A_ = -l E k = t, k + \g 2 + gNx 



0100) 2s = -1 A + = A_ = 1 E k = ^ + \g 2 - gNx 



1110) 2s = 1 A + = l A_ = E k = ?>^+\g 2 + gN. 



x 



0111) 2s = -1 A + = -l A_ = E k = 3Z k + ±g 2 -gN 



x 



1011) 2s = 1 A + = A_ = -l E k = ^ k + \g 2 + gNx 



1101) 2s = -1 A + = A_ = l E k = ^ k + \g 2 -gNx 



— one 4- dimensional subspace M k9 spanned by the vectors 

10000) = |1), |1001) = |2), |0110) = |3), |1111) = |4), 

and where 2s = A + = A_ = 0. If we choose for the basis of the M kg the following set of 
vectors: {|1) , — |2) , |3) , — |4)} and denote the projector on M k9 by P k9 , then the matrix 
form of the operator h k — C k is 



This leads to the secular equation 



2A fc \ 



/ 

2i k _ + g 2 

2£ k + g 2 

V 2A fc 44 / 



(32) 



(24 + g 2 - E) 2 [E 2 - 4£ k E - 4Aj] = 0, 



(33) 
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which yields the following eigenvectors and eigenvalues of the operator (1321) : 

|1001) 24 + g 2 

I0110) 24 + s 2 

u k _ 10000) + v k 11111) 2£ k -2E Gk 
u k |1111) - Vfc |0000) 2£ fc + 2E Gfc 



where 



2 V -Eg/c / ' - 2 V E Gk 



VI. THE STATISTICAL SUM AND ORDER PARAMETERS 

The operators h k and <t q act in different subspaces, so the statistical sum of the system 
equals 

z = n%n z - = ri Tr ex p i-p ih+m- 1 ^)] n^p ims^ 

fc>0 Q fc>0 a ^ 

= exp[-/3C ][2cosh(^a;| ( S'i)] iV J]; Trex P[-«- 

fc>0 

The eigenstructure of h k , found in Sec. [V] yields 

Z k = Tr exp {-(3h k ) = exp [-2/5 (A k r k + 4)] 8 exp (-§/fy 2 ) cosh(f3Q cosh(j3gy) 

" (35) 

+ 2 exp (-/V) (1 + cosh (2(3gy)) + 2 cosh {2j3E Gk ) + 2 , 
where y = Nx. Eq. (1281) can be rewritten in the form 

2r k = (3^^-\nZ k . (36) 

Using Eqs. (1271) . ( 1361) we get the equation for the order parameter 

A, = llAI-^G^-^-Fj (ftJS^fr.y) , (37) 
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where 

sinh(2/3£ G fe) 



F, (ft E Gkl 4, y) 



4e-/ ? 9 2 /2 cosh(/3gy) cosh(/?f fc ) + e"^ 2 [1 + C osh(2/fyy)] + cosh(2/3£ Gfe ) + 1 ' 

"(38) 

Passing from summation in Eq. (1371) over k to integration over the single-fermion energies 
£ one obtains for a sufficiently thin band S: 

s 



& = \G Q J ^F x %E G ,i,y)$i, (39) 



-<5 



where q denotes the density of states in S. Equation (|39|) is similar to the gap equation in 
BCS theory, however, the convexity properties of F^ft, E G , £, y), with respect to A, differ 
in general from those of F BCS ({3, E G ) = ta.nh(PE G /2y^. 

Given the spectrum of h Q , one easily computes the thermal average 

(%- " n k+ + a a ) ho . (40) 
From the commutation relations (l3Tj) one obtains 

(%- - %+ + °a) ho = ( n k- - U k+) hk + (Oha > ( 41 ) 

where 

/ N Tr (nfc_ - n k+ + - n_ k+ ) exp[-(3h k ] 

\ n k- ~ n k+) hk ~ 2 Tr exp J (42) 



and therefore 



Tr (n k _ - n fc+ + n_ k _ - n_ k+ ) exp[-(3h k ] (f3g) ^Tr exp[-(3h 



2 Tr exp [-/%] 2 Tr exp [-/%] 

The eigenstructure of yields 

Tr exp[-/?/g = exp [-2/3(A,r, + Q] [8 exp{-\(3g 2 ) cosh(/%) cosh(Pgy) 
+ 2 exp(-/3# 2 ) (1 + cosh{2(3gy)) + 2 cosh(2/3£ Gi J + 2] . 

Hence 

(n k _ - n k+ ) hk = [2exp(-±/3# 2 ) cosh(/3&) sinh(/fyy) + exp(-/V) sinh(2^/)] 
xsinh- 1 (2/5E G ,)F(/5,E G ,,4,y). 

Furthermore, 

Tr exp[— f3h a ] = 2 cosh^Pgyd^ 1 ), 



(43) 
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{a a ) ha = 2t^\i{(3gyd- 1 ). (44) 
Therefore, Eq. (TlBl takes the form 

y = F 2 (P,E Gk ,£ k ,y), (45) 

where 

2e-^ 2 / 2 cosh(g4) sinh(^y) + e^ 2 sinh(2^y) 

alP. Gfc> 2/J - ^ 4e -^ 2/2 cosh ( to ) cosh(^fc) + e-^ 2 [1 + cosh(2/fyy)] + cosh(2/3£ G J + 1 

+ gtanh^gyd" 1 ). 

(46) 

Equations fl39l) and (j45p constitute the set of two equations for A, y. The expression for the 
statistical sum (J341) and the equality 

resulting from Eq. (|36l . as well as Eqs. (|22l . (J23l) . lead to the following expression for the 
free energy 

T = min \ A 1 E gIFxW, E Gk) y) + y 2 + <? 2 + 2£ fe - /TMnP* 

n (47) 
-2^iV(|5|/3)- 1 ln[2cosh(/5^- 1 )] , 

fc>0 

where the minimum runs over all solutions of Eqs. (jSHD, (HB]) and 



Se"^ 2 / 2 cosh(/%) cosh(/fyy) + 2^ (l + cosh(2^y)) + 2 cosh(2/?£ Gfc ) + 2 



If G kk , is nonvanishing and constant only in a thin band near the Fermi surface: G kk , = 
Gx(k)x(k!) j C > 0, then the free energy density equals 

s 

f(h,P) = min e /" [lA^F^A^^^ + l^ + ^ + e-l^lnPld^ 

(«) 

- pd/T 1 y In [2 cosh^ycT 1 )] d£ + E (A = 0) + £<5 2 , 
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where d = NlS^ 1 is proportional to impurity concentration and E (A = 0) denotes the 
ground-state energy of free fermions. The last two terms are the contribution to the free 
energy density from one-fermion states, lying outside S. 

Further sections contain numerical analysis of Eqs. (159~j) . (|45[) and (|48p . The numerical 
solution of the set of Eqs. (1551) . (1431) has been carried under the assumption that one of the 
two inequalities 



at T = OK, whereas Eq. (1431) is satisfied by y = g in this limit. The above inequalities state 
weak coupling of impurities to conducting fermions. Thus at sufficiently low temperatures 
T, close to OK, A(T) is taken as the solution of Eq. (j49j) and substituted into Eq. (|43|) . which 
is then solved for one-fermion energies £ G S by exploiting the Newton-Raphson method. 
The resulting values of ?/(£) are used to obtain A(T + AT) from Eq. (1391) by deploying 
a Newton-Cotes quadrature as long as the result is self-consistent. The resulting value of 
A(T + AT) is exploited to achieve y(£) at temperature T + AT. This procedure is continued 
until T reaches the specified value. 

The computations in Sees. IVII1 IVIIII are performed under the assumption that /i = e F 
and dfi/dT = 0. Justification of this assumption is given in Sec. [DO 

VII. PHASE DIAGRAMS 

Equations (139|) . (1431) clearly possess the solution A = y = at all values of (3 > 0. At 
sufficiently large values of (3 one finds also other solutions, viz., {A ^ 0,y = 0}, {A = 
0,?/7^ 0}, {A 7^ 0,y 7^ 0}. The system's state is characterized, according to Eq. (1481) . 
by the solution which minimizes f(h,j3). It will be denoted by {A m ,y m }. Accordingly, we 
distinguish the following phases: 

— {A m = 0, y m = 0} corresponds to the paramagnetic phase P 



--g 2 + gy + Z<2E G 



or 



2gy -g 2 < 2E G 



holds at T = OK. Eq. (1371) then takes the limiting form 




(49) 
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— {A m 7^ 0,y m = 0} corresponds to the unperturbed superconducting phase S 

— { A m = 0, y m 7^ 0} describes the ferromagnetic phase F without bound quartets, in 
which the impurity spins tend to align opposite to those of conduction fermions 

— {A m 7^ 0,y m 7^ 0} describes the D phase in which superconductivity coexists with 
ferromagnetism 

We define the following temperatures corresponding to the respective phase transitions and 
characteristic points of A(T) plots: 

— T c , 1st or 2nd order transition S — > P, the order depending on g 

— Tpp, Curie temperature of 2nd order transition F — > P 

— Tsd, 2nd order transition D — > S 

— Tfd, 2nd order transition D — > F 

— Ti, T* are the end-points of the interval where A(T) is double- valued 



6,0x10 
5,0x1 4 - 

> 

^ 4,0x10" 4 

p 

< 

3,Dx10 
2,0x1 4 
I.OxlO" 4 



0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 
T [K] T 



7,0x10 
6,0x10^ 
5,0x1 J 
4,0x10" 
3,0x10" 
2,0x1 J 



- A (T) for y m > 

■ A(T) for y = 

■ A(T) for y > 



0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5 
T[K] 



(a) (b) 

FIG. 1: The gap parameter plotted as a function of reduced temperature T for Gq = 0.3,(5 = 
0.01 eV, d = 0.05 and: (a) g = 0.01 VeV, (b) g = 0.03 VeV. For sufficiently large value of g the 
phase transition S — * P is of the second order and T c = T* . 



The set of Eqs. (1391) . (H5l) has been solved numerically for different values of g and 
impurity density d. FIG. [1] shows how variation of g affects A(T). For solution with 
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(a) 




(b) 



I Pure SC phase 

I Perturbed SC phase 



■ Impurity phase Gapless region 

l Paramagnetic phase 






T[K] 



T[K] 



(c) 




(d) 




T[K] 



T[K] 



(e) 




(f) 




T[K] 



T[K] 



FIG. 2: Phase diagrams of the system for the following values of parameters: Gg = 0.3, 5 = 0.01 eV 
and: (a) d = 0.001, (b) d = 0.001 (detailed version of (a) for small values of g 2 /d) , (c) d = 0.01, 
(d) d = 0.01 (detailed version of (c)), (e) d = 0.1, (f) d = 0.5. The gapless region denotes the zone, 
where the smallest excitation energies from ground state disappear. 

A ^ Oj = 0, \T* — T c \ declines as g is increased and finally vanishes. The temperature 
range, where A(T) is nonvanishing declines for the solution {A ^ 0,y ^ 0} with increasing 
value of g and is smaller than for the solution {A ^ 0,y = 0}. 
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The phase diagrams for fixed values of Gg = 0.3, 5 = 0.01 eV are depicted in FIG. El 
These diagrams show the decline of the D phase with increasing impurity concentration d. 
The S phase and the T c temperature, on the other hand prove to be almost insensitive to 
variations of d, but one notices increase of T c with g 2 N~ l oc g 2 d~ l at higher d. In this 
respect the system's behaviour is rather unusual and may be explained as a consequence of 
the weakness of V and breakdown of Eq. (|22|) at larger values of d. However, it is worth 
noting at this point that solid solutions of magnetic metals in zirconium and titanium have 
higher superconducting transition temperatures than pure zirconium and titanium^. 

A subregion NG with gapless superconductivity and A m > (the light gray strip in 
FIG. [2]) is present for sufficiently large values of g 2 d~ 1 . The appearance of NG subregion is 
due to the negative terms — gy, —2gy, which contribute to the spectrum of h^. 

The perturbative effect of impurities increases with g 2 d~ 1 and causes the complete sup- 
pression of superconductivity at low temperatures, when g 2 d~ x is sufficiently large. At 
sufficiently low temperatures the mixed phase is present between the S and F phase, where 
superconductivity coexists with ferromagnetism of impurities. 

The temperature Tpp decreases with impurity concentration d and depends almost lin- 
early on g 2 d~ x . One also observes decrease of Tpp, with g 2 d~ l and increase of T$d with 
g 2 d~ 1 . The temperatures T S d, T c , T fd , T pf satisfy the following inequalities: 

Tsd < T c < Tpp for small values of g 2 d~ 1 

Tfd < Tpp for large values of g 2 d~ 



VIII. SPECIFIC HEAT UNDER VARYING IMPURITY CONCENTRATION 



In terms of the energy density u = ^M- the system's specific heat expresses as 

du du dy du dA 



+ 



+ 



dT dy dT dA dT 



From Eq. (1481) one obtains 



(50) 



6 

u (P, E g , £,y)=ej [| A 2 E G X (F, (P, E g , £,y) + /3 



dF\ 



(51) 
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The specific heat corresponding to a definite phase is denoted by the appropriate index, e.g. 
cp denotes the specific heat corresponding to paramagnetic phase P ({A m = 0, y m = 0}). 



0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 
T [K] T 




-c n (T)/c p (T) ford = 0.5 
-c D (T)/c p (T) ford = 0.3 
-c n (T)/c p (T) ford = 0.1 
-C E (T)/c p (T) 
-c„(T)/c p (T) 



4,0 4,5 5,0 



(a) (b) 

FIG. 3: The specific heat for different values of impurity concentration. (Gg = 0.3, 5 = 0.01 eV 
and (a) g^d' 1 = 6.25 x 10" 5 eV, (b) g 2 d~ l = 2.50 x 10~ 4 eV) 

The specific heat ratios c(T) /cp(T r ) under varying impurity concentration d are depicted 



in FIG. |3j Similarly as in Ref. [14| no specific heat jump Acsp = c$ — cp at T c for sufficiently 
small values of magnetic coupling constant g is observed and very weak dependence of cs 
and cp on d has been found. FIG. [3] also shows how the S phase is replaced by the D phase 
at low temperatures as d increases and how the jump of c at Tsd increases with d. 



IX. CHEMICAL POTENTIAL 



The properties of a superconductor can be determined from the solution of the gap 
equation, which is supplemented by the following equation for the chemical potential \x 



(52) 



n denoting the average number of fermions in the system and po is the trial density matrix. 
For a superconductor with 4-fermion attraction perturbed by magnetic impurities po equals 

exp(-/3/i ) 



Po 



Tr exp(— f3h ) 
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and Eq. (152]) takes the form 

Y Tr ( n k*po) = \ Y Tr + n -k+ + + n ^-)Po] = -li 3 ' 1 Y 



k>0 



fc>0 



IE 

fe>0 



G'A' 



Eaki 2/) — 2F 3 (/5, E Gk} £ k , y) 



dZ k 
<% 

- n. 



(53) 



Thus, in the thermodynamic limit 



1 - ^-F t ({3, E G) £, y) - 2F 3 (P, E G , £, y) 

&G 



d£ = d, 



(54) 



where 

F 3 (f3,E G ,£ k ,y) 



sinh(/^) cosh(Pgy) 



4e"^ 2 /2 C osh(/%y) cosh(/3£ fe ) + e"^ 2 [1 + cosh(2/fyy)] + cosh(2/3£, 



Gk) 



(55) 



and (i denotes the density of fermions in the free electron gas model, viz., d = ^£pQp, Qf 
mp F {2n 2 K 2 )~ l . 

The equation for \i at T — results from Eq. ( |54l) in the limit T — > 0, viz., 



<?(£) 



-Mo) 



E, 



G 



d£ = d. 



(56) 



Substituting E G = (£ 2 + A 2 (0)) 1 / 2 , g(£) = Q F x(0 for £ G [—5, 5] and exploiting the equality 



q(£) = QfvZ + / i (0)/v / ^ one obtains the equation for //(0) [cf. Ref. [l4j 

|^ 1/2 (/i(0) - 5) 3 / 2 + 25 = \e F . 



(57) 



For 5 = 0.01 eV, e F = 1 eV, Eq. flSTD has the solution /x(0) = 0.9999748 eV. For 5 = 0.01 eV, 
e F = 1 eV, g = 0.01 v 7 ^, d = 0.05, A(T X ) = y(Ti) = 0, Eq. ([54]) has the solution /i(Ti) = 
0.9999754 eV, which is the same as /x(T*) for A(T*) ^ 0, y(T*) ^ 0. It follows that the 
relative difference |/j,(0) — /x(Ti)|//x(Ti) is of order 10~ 7 and |/u(Ti) — £p| is of order 10~ 5 eV. 
The very weak variation of fi(T) does not affect the temperature variation of A(T), y, free 
energy and specific heat within the accuracy applied in numerical calculations. 



X. CONCLUSIONS 



We have shown that the thermodynamics of a superconductor with a quartet binding 
potential, perturbed by a reduced s-d interaction is solvable if the impurity density is suffi- 
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ciently small. The essential properties of a superconductor perturbed by magnetic impuri- 
ties have been demonstrated: decline of the temperature range, where the gap parameter is 
nonvanishing, with increasing coupling between conduction fermions and impurities, strong 
dependence of specific heat jump on the density of impurities at the transition to the D 
phase with coexisting superconductivity and ferromagnetism and existence of gapless super- 
conductivity. Due to the assumed weakness of the s-d coupling the value of A m (0) and T c 
are not diminished by the impurities. These investigations will be extended to include the 
effect of a general s-d exchange interaction V s _ d . 
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